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Abstract 

We apply different integrability analysis procedures to a reduced (spatially homogeneous) mechanical 
system derived from an off-shell non-minimally coupled N=2 Maxwell-Chern-Simons-Higgs model that 
presents BPS topological vortex excitations, numerically obtained with an ansatz adopted in a special 
- critical coupling - parametric regime. As a counterpart of the regularity associated to the static 
soliton-like solution, we investigate the possibility of chaotic dynamics in the evolution of the spatially 
homogeneous reduced system, descendant from the full N=2 model under consideration. The originally 
rich content of symmetries and interactions, N=2 susy and non-minimal coupling, singles out the proposed 
model as an interesting framework for the investigation of the role played by (super-) symmetries and 
parametric domains in the triggering/control of chaotic behavior in gauge systems. 

After writing down effective Lagrangian and Hamiltonian functions, and establishing the correspond- 
ing canonical Hamilton equations, we apply global integrability Noether point symmetries and Painleve 
property criteria to both the general and the critical coupling regimes. As a non-integrable character 
is detected by the pair of analytical criteria applied, we perform suitable numerical simulations, as we 
seek for chaotic patterns in the system evolution. Finally, we present some Comments on the results and 
perspectives for further investigations and forthcoming communications. 



1 Introduction and Motivation 

The dynamical properties of gauge systems define a focus of great interest, as one can realize from the 
remarkable effort that has recently been driven to the analysis of the stability of gauge field configurations 
[1]. As a cornerstone of theoretical physics, gauge theories have been intensively investigated, and the actual 
result is a comprehensive, but non-exhaustive, picture of such a successful theoretical framework. One 
(promising) aspect of such systems is the room for coherent soliton solutions, that may play a crucial role in 
the understanding of physical phenomena like, for instance, the quark confinement. On the other hand, the 
search for chaotic regime windows in gauge theories seems to be as much important as the former analysis, 
defining a counterpart approach that may lead to answers to key-problems, for instance, and again, the 
confinement phenomenon [1, 2]. In the eighties, a method for investigating chaos in field theories has been 
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developed and applied to Yang-Mills systems [2, 3]. The main idea is to reduce the model to its mechanical 
limit, by considering spatially homogeneous field configurations. The discussion of chaotic evolution in this 
restricted regime is conjectured to be sufficient to ensure chaotic behavior for the full field theory [2, 4]. 

In this context, a relevant focus of attention is defined by the possibility of establishing a systematic 
interdependence recipe for the relation between gauge symmetries and the control of chaotic dynamics. 
Among gauge symmetries, one should eventually care about super symmetric systems, considered either as 
manifestations of a fundamental symmetry or as enriched models conceived to be a tool to better describe 
physical situations. 

Planar (2+1) Maxwell-Chern-Simons-Higgs (MCSH) theories, as candidates for an effective description of 
high-T c superconductors phenomena, have recently been chosen as target models to the aim of order-to-chaos 
transition studies. Bambah et al. [5] have considered both the (proven to be) integrable minimally coupled 
Chern-Simons-Higgs (CSH) model and its higher momenta natural extension, namely, the minimally coupled 
Maxwell-Chern-Simons-Higgs system. The latter failed when submitted to an intcgrability criterium, the 
Painleve test, leaving room for a chaotic regime that happened to be confirmed by numerical Lyapunov 
exponents and phase plots analysis. Recently, Escalona et al. [6] have performed a similar work upon 
a MCSH system endowed with both minimal and non-minimal couplings in the interaction sector. The 
non-minimal coupling stands for a Pauli-type term describing a field-strength/matter-current interaction, 
admitted in (2+l)-D regardless of the spin of the matter field [7]. Moreover, if quantum extension is a claim, 
such a non-minimal coupling should be considered from the start [8] . In the work of Rcf. [6] , the CSH system 
is argued to be still integrable, while the non-minimal MCSH exhibits "alternating windows of order and 
chaos", as the non-minimal coupling constant g is varied, the other parameters defining a set of constant 
inputs. The model they adopt is the bosonic projection of an already established N=l- super symmetric 
system. As a matter of fact, a non-minimally interacting MCSH system had formerly deserved an extension 
endowed with on-shell N=2-susy[9]. As far as soliton solutions are a subject of interest, the N=2 extension 
defines the proper framework, allowing for the self-dual regime [10, 11]. On this token, Antillon, Escalona 
et al. have found, in a previous work [12], while working upon such an N = 2 extendable model, a self-dual 
static non-topological vortex solution, motivating their search for spatially homogeneous chaotic dynamics 
as an interesting counterpart for their former discovery. Nevertheless, even if one assumes the validity of the 
conjecture relating the mechanical limit to the full theory, a problem arises if the counterpart character is 
claimed to be rigorous: the vortex has been found in an N=2-susy framework, while the varying g procedure 
adopted in [6] necessarily moves the system out of the N=2-susy-bosonic projection situation. As clearly 
assumed in [9], a critical coupling, namely, g = —e/n, has to be verified to ensure on-shell lSl—2-susy, where e 
is the minimal coupling constant and k is the Chern-Simons mass parameter. Moreover, the scalar potential 
is forbidden to be anything but the non-topological mass-like (j> 2 term. So, varying g while keeping e and k 
constant, and adopting V = \{(f) 2 — v 2 ) render their model a sector of, at most, an N=l system. 

Alternatively, another planar N=2 non-minimal MCSH model has been recently proposed [13], defining 
a richer spectrum that presents both non-topological and topological self-dual static vortex solutions [14], 
numerically obtained after the adoption of the critical coupling relation. Such a system exhibits off-shell- 
realized N=2-susy, and is obtained from a N=l-D=4 ansatz, after dimensional reduction and a suitable N=2- 
covariant superfield identification. Two important differences arise, if one settles a comparison encompassing 
both non-minimal MCSH models: in the N=2-off-shell case, an "additional" 2 neutral scalar field takes 
place; also, in the N=2-off-shell case, no relation between the coupling constants and parameters is required 
to ensure N=2-susy (though the vortex excitations have so far been shown to prevail in the particular 
g = —e/n regime). In other words, if the model of Ref.[13] is considered, the freely varying g strategy and a 
topologically non-trivial scalar potential happen to be compatible with N=2-susy. 

2 A common, improper, terminology. The "minimum" content forbids interesting excitations, like topological vortices. 
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Motivated by these interesting features, we carry out the analysis of the reduced - mechanical - version of 
the bosonic-sector Lagrangian extracted from the Rcf.[14]. In the next section, we present the theory, the field 
equations and their spatially homogeneous counterpart, the one-dimensional effective Lagrangian and the 
associated conjugate momenta. As we detect an additional (besides the Hamiltonian) constant of motion, a 
convenient reparametrization is implemented, and the corresponding Hamiltonian system is displayed, ending 
up with the canonical Hamilton equations. In Section 3, we move back to a (general regime) second-order 
formulation, upon which we start our analysis of the integrability issue, adopting two alternative analytical 
criteria - the Noether point symmetries approach [15] and the Painleve test procedure [17, 18]. The former 
strategy leads to a set of ten coupled partial differential equations which seem to possess a closed form 
solution only in the minimal coupling regime, g = 0. In the Painleve test context, it is shown that no set of 
negative integer dominant exponents can be found, as the equation for the gauge field is considered, spoiling 
the corresponding test algorithm for the verification of a possible strong Painleve property associated to the 
system we propose. The role played by the gauge field equation suggests that a change in the gauge sector 
dynamics might render a different situation for the integrability analysis. So motivated, in Section 4 we 
adopt the critical coupling relation regime, arriving at first-order equations for the gauge degrees of freedom. 
We present both the associated effective Lagrangian and Hamiltonian functions, as well as the corresponding 
iterated second-order field equations. These results enable us to re-address the Painleve test, which indicates 
that the critical coupling regime presents an even worse feature concerning the presumed strict negativeness 
of the dominant exponents. Revisiting the Noether point symmetries approach also gives no clues on possible 
integrable setups. In Section 5, an analysis of chaos is performed with physically acceptable values of the 
parameters; in the regimes characterised by g ^ 0, for both non-critical and critical values of g , the model 
becomes more stable than in the g = situation. Finally, we present our Concluding Comments and we 
discuss the failure of the system in obeying the strong Painleve property, a fact that may lead to future 
consideration of other integrability tests. We also try to explain why g ^ yields a more stable behavior in 
spite of the inherent non-linearity brought about by a non-vanishing g. 

2 Describing the Model 

We start [14] with: 

-l(d,M)(d^\ 2 ) + ^A^-U, (1) 

where 

and V M = (d/j, — ieA^ — igF^)4>. G is defined as G = 1 — g 2 \4>\ 2 - The field equations for the full theory 
read: 

d^p + mpp = -j p -Vv^ 

where = f {<j)*V ^(j) - 0(V M 0)*), and also 

d^M - l da r\tf + e 2 (l^| 2 ^ 2 + (V e )M + 2 g M|0| 2 )(2./ e + 2 g |0| 2 ) + e2 ^ 2M = Q 
2 4G- 
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and -(V a V Q 0)*-|^(a Q a Q M) + 



|0| 2 - v 2 + (2/s/e)M + 2 9 |^)| 2 M 



8G 2 



- + 



+ - 



|</>| 2 -v 2 + {2n/e)M + 2g\<j)\ 2 M (1 + 2<?M) 



J- A* 



AG 



+ ■ 



M 2 



= 



Adopting the gauge choice A = and imposing the spatial homogeneity, namely, <9j(V field) = 0, the 
phase of the scalar field becomes a variable with vanishing time-derivative, and one can eliminate it without 
loss of generality. So, one ends up with a real scalar field and the following set of reduced equations of 
motion: 



d_ 
~dt 

d_ 

dt 



GAi 



-2e#— [<j>A 2 ] - e 2 <t> 2 A 1 - kA 2 



+2eg(j)— [4>Ai] - e 2 (t) 2 A 2 + kA x 



(2) 
(3) 



\ , = 1(62) _ e 2 6 2 M _ e 2 (0 2 - v 2 + (2«/e)M + 2gM^ 2 )(2^/e + 2g<j> 2 ) 
2 4G 



-(eAi + ,gA 2 ) - (eA 2 - + g4>M - e 2 <jjM 2 

e 2 g 2 4>{4> 2 - v 2 + (2/c/e)M + 2gM<j) 2 ) 2 
4G 2 

_ i; 2 + (2«/e)M + 2gM0 2 )(l + 2gM) 
2G 

The effective Lagrangian that generates these equations of motion is: 



r,2±(±2 



L = 



G 



(A,) 2 + (A 2 f] - I (A,A 2 A 2 A,) e ^f (A, 2 + A 2 2 ) + 



(My 



g4>(j)M 



(<j>f e 2 2 - v 2 + (2 K /e)M + 2gM(f> 2 ) 2 e 2 (f> 2 M 2 
~~2 8G 2 



where G =1 - g 2 <j) 2 and Q = k + 2eg4> 2 . 

The canonically conjugate momenta, defined, as usually, by p — have the expressions: 



7Tl 



dL 
dL 

TT 2 = — 

dA 2 
dL 

P<P = — y 



M 



dL 
~d~M 



GA 1 + |a 2 , 
Gi 2 - I Ax , 
4> - 9<t>M , 
M - 5 # . 



Before we proceed to the canonical Hamilton equations, let us notice that the quantity 

I = A 2 TTl - AlTT 2 



(4) 
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is a constant of motion. This can easily be checked by multiplying and combining the "gauge fields" equations 
according to ^(2) — A\(3). Motivated by this fact we reparametrize the gauge sector adopting polar 
coordinates, instead of Cartesian ones. We have: A\ = Acos(, A 2 = Asen(, and the "new" set of variables 
is (A, (, <j), M). The Lagrangian now reads: 



(i) 2 + ^( C ) 2 ]-|(^) + ^-^M + 



(M) 



e 2 <j> 2 A 2 _ e 2 (0 2 -v 2 + (2«/e)M + 2gM^ 2 
2 8G 



3 Vm 2 



(5) 



yielding the same expressions for and Pm, and defining Va = ^j l = GA , p^ = ^ = GA 2 ( — ^-A 2 . One 
can easily check that = —I, resulting p\ = 0. The Hamiltonian reads: 



H C AN. = ^ 



M 



+ 



1 

+ 2G 



A 2 + 



(Q/2) 2 + e 2 G(b 2 

The canonical Hamilton equations result: 

PA 



+ ^(<f- v2 + (2«/e)M + 2g4> 2 MY 



A = 



G 



PA = 



G 



pT 

A 3 



- ( K (Q/2f + e 2 G^A 



1 

G 



P£,Q 
A 2 2 



PC = 0; 



(6) 



— [p + g<f>P M ] ; 



P<t> = - 



G 2 



2 , P( , 2 , n 2 
PA + J^2+P4 +PM 



■ + 2e)p c + g(l + g 2 <f) 2 )p 4> P M 



iM 2 



-v 2 + (2n/e)M + 2gM<j> 2 )' 



Hl , !l + 2ey^+e 2 G 2 <f>M 2 + ^ 

' ' " ' ' i> 2 - v 2 + (2«/e)M + 2gM<t) 2 ) (1 + 2gM) 



+ 



M = q[P m + 9<to>4] ; 



Pm = ~g 



e 2 G<f> 2 M + e l(4> 2 -v 2 + {2n/e)M + 2gM<f) + 2gcf 



3 Integrability Analysis: General Case 



We present two Lagrangian analytical criteria to address the issue of integrability: Noether point symmetries, 
better suited for establishing the constants of motion, and the Painleve test, meant to check for an overall 
property (the dependent variables being meromorphic for movable singularities on the complex time plane) 
that indicates integrability. 



G 



3.1 Noether point symmetries 

An important issue regarding a Lagrangian system concerns its Noether point symmetries, linking symmetries 
of the action functional to conserved quantities. Here we address the question of the existence of Noether 
point symmetries of our system following the method shown in reference [15]. 
We seek for infinitesimal point transformations of the form 



A 


= A + eijA, 


(7) 


c 


= ( + ?vc> 


(8) 


4> 


= 4> + £V<t" 


(9) 


M 


= M + erjM , 


(10) 


i 


= t + ST , 


(11) 



for rjA , ilc > , iJ<P an d t functions of the fields and the time, and e an infinitesimal parameter. These in- 
finitesimal transformations leave the action functional invariant up to the addition of an irrelevant numerical 
constant if and only if the following Noether symmetry condition [15] is satisfied, 

dL dL dL dL dL . -,dL 

T ~dt + VA dA +T1 <dt +T] *dj +riM dM + {r,A ~ ~ TA) dA + 

+ tL = F. (12) 

for F a function of the fields and the time. If such a function can be found, there is a Noether point symmetry 
and an associated Noether invariant I given by 

dL dL dL dL -dL -dL -dL ■ dL 

I = " A U + ^ + ^ + m M- T{A M + c ^ + ^ + M M- L) - F - (13) 

In the Noether symmetry condition, the time derivatives are to be understood as total derivatives, e.g., 

8t ■ dr ■ dr ■ dr • dr 

T= dA A+ dc c+ d-/ + dM M+ m- (14) 

Inserting in the symmetry condition the Lagrangian given by equation (5), we obtain that a cubic 
polynomial in the velocities must vanish. The coefficients of equal powers of velocities vanishing, we obtain 
a coupled set of linear partial differential equations determining both the symmetries and the Noether 
invariants. The cubic terms yields simply 

r = r(t) 1 (15) 

that is, the transformed independent variable is a function of time only. The equations associated to quadratic 
terms, however, are a set of ten coupled partial differential equations which, apparently, do possess a closed 
form solution only in the minimal coupling case g = 0. Restricting the treatment to this almost trivial 
case, and proceeding to the first and zeroth order terms, we just found time-translation and ( translation 
symmetries. These symmetries are associated, according to (13), to the energy and conservation laws. 
These are almost obvious results, showing that the nonlincarity and coupling in the potential gives no much 
space for the existence of conservation laws of the system, even in the g = case. This is a signature of 
non-integrability. However, other methods for investigating conservation laws of the system like Lie point 
symmetries [16] were not used. 
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3.2 Painleve test 

We now go back to a second-order configuration space formalism in order to settle the framework for the 
application of Painleve test [17, 18]. The equations read 

G 2 A = 2g 2 G^A+^-(^) 2 A-e(Kg + e)^A, (16) 
where C stands for the constant , 

G 3 4> = 9 2 G 2 ^ 2 -g 2 G 2 <pA 2 - g ^l-gC(Kg + 2e)(j>-- A {Kg + 2ef^A 2 + 
+ E -j(<P 2 ~ v 2 ) [eg 2 v 2 -2( K g + e)+ g 2 {2ng - e)cf + 2e fl V] 

- (ng + e) [k - ev 2 g + (3eg - ng 2 )cj> 2 - 2eg 3 cf> i + {ng + e)M] <j>M (17) 

and 

G 3 M = gG 2 4> 2 - g 3 G 2 cj ) 2 A 2 + - 9 -(k g + 2ef A 2 - g {k g + ef M 2 cf + 

+ ^(</> 2 - v 2 )[-2n + g (-4e + eg 2 v 2 + 2k 5 )</> 2 + 3eg 3 4 ] 

+ [-K 2 + (n 2 g 2 + eng 3 v 2 - e 2 - 'Seng + e 2 g 2 v 2 )<j) 2 + eg 2 {2ng - e)</> 4 + e 2 g 4 <j> 6 ]M . (18) 

Assuming time to be a complex variable, the first step of the Painleve test is concerned with the lead- 
ing singularity behavior. One supposes the leading terms to be of the general form A <~ a(t — to) a ,(j) ~ 
b{t — to) 13 ,M ~ c(t — to) 7 , where a,/3, 7 < 0. Such an assumption turns the last three equations into the 
following asymptotic (t — > to) relations: 

3 4 a6 4 a [a - 1 + 2/3] r Q+4/3 - 2 ~ ; 

,g 6 6 7 /3 [2/3 - 1] r 7 "" 2 ~ ff 6 6 5 a 2 a 2 r 2a+ ^- 2 + 
~2eg 3 ( K g + e)b 5 CT 5 ^ + {Kg + efbc 2 T 0+2 ^ ; 

C7(7 _ 1)r 4/3+ 7 -2 _ 9 1?/3(2I3-1)tW- 2 , 

where r = t — to- 
Starting from the last equation, one gets 7 = 2/3 and c = gb 2 /2. Inserting 7 = 2/3 in the second equation 
reduces the balancing to the first two terms, so leading to a = (3 and a 2 = (2/3 — l)b 2 / (3. But the first 
equation shows the impossibility of having a, (5 < 0, as a + 2/3 — 1 = is required, spoiling the Painleve test 
procedure. 

Another possibility would be to set a = in the first equation, leaving it behind as an identity. One 
could then drop the second term (first on right-hand side) of the second equation, and the balancing of the 
remaining three terms would lead to an interesting set of negative values for 7 and f3: 7 = —4, (3 = — 1, 
provided that the following relation holds: 



{ K g + efc 2 - 2{ng + e)eg 3 b i c + 3g 6 b 6 = 
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Still, one has to deal with a zero "dominant" exponent, which spoils the Painleve test. 

As one faces a problem with the gauge sector dynamics, the adoption of the critical coupling relation 
(projecting the system onto the regime that hosts the already established vortex excitations) may serve as 
a valuable tool of investigation. In fact, imposing g = —e/n leads to first-order equations for the gauge field 
[14]. 



4 Critical coupling regime 

If g = — e/K, one gets: 

kF v \Jv , 

and the reduction to spatially homogeneous configurations yields 

KGEijAj = -e 2 A l( j) 2 , (19) 
where £ i2 = +1 = —£21- From this set of equations one can arrive at 

and, as far as G > (a condition inherited from the original N=2-susy framework), this implies that A\ + A 2 , 
is a constant of motion (thus reproducing the "pure" minimally coupled Chern-Simons-Higgs situation). 
Adopting polar coordinates, A\ = Ccos(, A 2 — CsenQ (A 2 + A\ = C 2 ), and manipulating the set (19), 
one finds ( = —e 2 <p 2 /nG. Following the same route chosen in the general (non-critical) case, we seek for 
the effective Lagrangian and Hamiltonian, settle the canonical equations of motion and, as we aim at the 
Painleve test for integrability, iterate them to get second-order coupled differential equations. One can easily 
verify that the following Lagrangian and Hamiltonian functions 



L = 



PC 
M 2 



( + 



e 2 <j) 2 
kG 



- gWM - 



2 

s 2 (0 2 



+ 



2 2k 2 G 
'J 2 + (2K/e)M + 2g<j> 2 My 
8G 



e 2 M 



2^,2 



He AN. 



2G 



A 2 



+ QPC + P4> 2 + Pm 2 + 2g4>p<pPM 



+ 



+ 



2G 



(Q/2) 2 + e 2 Gtf 



A 2 + 



v 2 + {2n/e)M + 2g(j) 2 My 



where g = —e/K, G = 1 — (e 2 / ' n 2 )cf) 2 , and 
iterated second-order system turns out to be 



2 8G y 

= —kC 2 /2 3 , leads to the proper set of field equations. The 



-G A 



2/^2, 



-g G 



• 2 , g 2 C 2 <^ 



A 2 



+ [eg(C + egv^/A)} <f> + | [-2eg 2 v 2 ] cf> 3 + e 2 g 



1 + 



p 2 n 4 



3 In fact, such a relation between the conserved quantities and A 2 (= C 2 ) must be imposed to ensure first-order £ 
-e 2 4> 2 /kG equation of motion. 
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G 2 M = (g - 2g 3 <f> 2 + g 5 ^)4> 2 + (g<f> - 2g 3 ^ 3 + g 5 <f> 5 )4> + 2e 2 cj ) 2 M 

2„2,4», f 2,, f Ke ,2 2 2/-, , 4 , e V0 6 , ^V 2 



-e 2 g 2 ^M - k M - f<f - e 2 g 2 (l + ^-)^ + + 

Again, for the Painleve test, the asymptotic relations are found: <j> : (3 = or (3 = 1/2. If one takes the 
[3 = case, one is left with two problematic outputs, as the equation for M is considered: either 7 = (3 = 0, 
or (3 = 0, 6 2 = l/g 2 , 7 undetermined. So the signature of lack of strong Painleve property remains. 



5 Analysis of chaos 



Since the results of the analytical approaches suggest that the system may not be integrable, we now turn 
to a numerical study to verify if such a non-integrability feature is presented in a chaotic form. 



5.1 SALI method 

The most well-known method used to detect whether a system is chaotic or not is the maximal Lyapunov 
Characteristic Exponent (LCE), a\. If g\ > the flow is chaotic. The G\ is computed [19, 20] from 

L t = ^ In |~=^jyj i performing the limit (20) 



a x = lim L t , (21) 

t — *oo 

where w(0), w(t) are deviation vectors and the time evolution of w is given by solving the equations of 
motion and associated variational equations. 

Since these vectors tend to acquire an exponential growth in short time intervals, many calculations of 
Lti, as w(t) evolves for a short time t\, are carried out after each w(t) is normalizated. With this procedure, 
the mean value of Lt x is computatcd as 




i=l 



For Hamiltonian sytems, this computation becomes very lengthy with poor convergence, and this long 
procedure may point to a false chaos diagnostic. 

We have chosen to adopt the method developed by Skokos, Antonopoulos, Boutis and Vrahatis, the 
so-called Smaller Alignment Index, SALI, for brevity [21, 22]. The reason for this choice is that the SALI 
method is computationally faster and less unstable than the Lyapunov exponent analysis, improving the 
adequacy of the former for the system we investigate. The SALI is a indicator of chaos that tends to zero 
for chaotic orbits, while it exhibits small fluctuations around non-zero values for ordered ones. So, the SALI 
is defined as: 



SALI(t) = min 



Wi(t) w 2 (t) 



+ 



wi (t) w 2 (t) 



\\m(t)\\ \\Mt)\ 



(22) 



\\m(t)\\ \\Mt)\\ 

where W\(t) and w^t) are the evolutions of two deviations vectors with different initial conditions, || • || is 
the Euclidean norm and t is the time . 

The authors of SALI method showed that SALI can be approximated by means of the difference of the 
two largest Lyapunov characteristic exponents u\ and er 2 - 
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The main advantage of the SALI in chaotic regions is that it uses two deviation vectors and exploits 
at every step the convergence from all previous steps. The SALI value tends to zero for chaotic flows at 
a rate which is a function of the difference of the two largest Lyapunov characteristic exponents ui, a<i as 
SALI oc e - ( <Tl_cr2 ) t . As usually done in numerical computations, we need to define a threshold so that a 
computed number be considered zero. In most of the cases, the selected value is < 10~ 5 . Like in the case 
of the Lyapunov exponent, it also happens that, in the SALI method, this is the criterion we shall use to 
distinguish between order and chaos. 



5.2 Equations of motion and requirements 

The integration of the system was performed by means of a Gear algorithm, in a variable step mode, starting 
from the minimal step size h = 0.0001, and eventually getting reduced in order to preserve the value of the 
Hamiltonian and pq , known to be constants of motion. Another constraint maintained along the integration 
was that G (G = 1 — g 2 <j> 2 ) should be greater than zero. The following first order equations of motion were 
used in the numerical integration: 



A - V -± 
G 



1 

PA= G 



A 2 



^-([Q/2f+e 2 G^)A 
PC = 0; 



— [p + g<j>P M ] ; 



PA 2 + -jt + P<p 2 + Pm 2 



+ {ng + 2e) 



A 2 



+ g<j>{Kg + 2e)p c + g(l + g 2 ^) V4> P M 
+ e 2 G 2 4>M 2 + (0 2 - v 2 + {2n/e)M + 2gM<f> 2 ) 2 



+ 



e 2 G(f> 



(<t> 2 -v 2 + (2n/e)M + 2gM<f> 2 ) (1 + 2 5 M)| ; 



Pm — —- 



G 



M = - [P M + g<pP4>] ; 
e 2 G<p 2 M +^(4> 2 -v 2 + (2«/e)M + 2gM<p 2 ) + 2# 2 



For the sake of simplicity, we have adopted the following notation: 

C = Qi, A = q 2 , 4> = q 3 , M = q 4 , p ( =pi,p A = Pi, P4, = P3, Pm = Pi- 

For each set of parametric inputs the numerical integration was performed and SALI method used after 
transient damped. In the following, we display representative samples of our findings. 
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5.3 Case with g=0 



Since our model comes from a supersymmetric version of Maxwell-Chern-Simons-Higgs with non-minimal 
coupling, it is not clear whether we shall recover dynamical properties similar to the ones observed in other 
studies [6, 5] when g = is adopted, and the model is so reduced to the minimal coupling case. The two 
main properties found in these studies were the existence of chaos in the presence of Maxwell term and the 
asymptotic evolution - order versus chaos - sensibility to initial conditions. To verify these properties in our 
model, we have chosen the initial conditions defining a fixed point of the system, and then we varied qs from 
to 2 with parameters set as e = 2, k = 2, v = 2, and as mentioned above, g = 0. The initial conditions are 
gi (0) = 1, q 2 (0) = 1, qi(0) = 2, Pl (0) = 0, p 2 (0) = -1, p 3 (0) = and p 4 (0) = 0. 

The results are presented in a sequence of phase diagrams. In Figure 1, we display p4 versus q±. We 
remind that we could have used any pair of variables to perturb the initial fixed point and then present the 
outcome in these graphs, but our choice was guided by the fact that both variables q% (<fi) and q4 (M) are 
present in the Higgs-type potential of our model. After this, in the Figure 2, we show a graph of SALI as a 
function of q 3 . 




a)g 3 (0) 



= 0.1, b)« 3 (0) = 0.3, c) g 3 (0) - 0.7, d)<a(0) - 1.5, e) g 3 (0) = 1.75 and f) q 3 (0) = 2.0. 

Figure 1 
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Figure 2 



We can see from Figure 2, with the help of the break in the vertical axis, that for q 3 > 0.55 the behavior 
of the system becomes chaotic with SALI assuming values between 10 -89 and 10~ 33 . 



5.4 Case with g ^ but outside critical coupling regime 



Now, we use g ^ 0, but outside critical coupling regime, to check whether the inclusion of this kind of coupling 
may turn some configuration dynamics into chaotic, with initial conditions that, in the case g = 0,lead to 
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regularity. To do that, we fix the parameters and the initial conditions as e = 2, k = 2, v = 2 and ?i(0) = 1, 
q 2 (0) = 1, q 3 {0) = 0.25, q 4 (0) = 2, pi(0) = 0, p 2 (0) = -1, p 3 (0) = and p 4 (0) = and we analyze the model 
with the following values of g: 0.1, 0.7, 1.5 and 2.5. 

For all these cases, the behavior remained the same, indicating that the variation of g does not change 
the behavior of the system from regular to chaotic, as we can see in Figure 3. 




I — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ Ml i . , ■ , . i ■ i . i i — i— 

4fi 4A A3 0,0 0,2 0* 0,» -0* -0,4 -0^ 0.0 0,2 OA bfi 

P 4 



a) 5 =0.1,b) 5=0.7,c) 0=1.5 and d) 5=2.5. 
Figure 3 



5.5 Case with g in the critical coupling regime 

Now, we explore the critical coupling regime where g = — — . We fix the parameters k — 2, v — 2 and vary 

k 

e , keeping the initial conditions as a perturbation case of the fixed point; but, in this case, we shall have a 
different set of initial conditions for each e, since the general expression for fixed point element p4(0) depends 
on e. With this in mind, we keep the same values for qi(0) = 1, 92(0) = 1, Pi(0) = 0, ^2(0) = — 1, ^3(0) = 0, 
P4(0) = and set 93(0) = 0.7, a value that makes the system chaotic in the case g = and 54(0) = e. In this 
case, (73(0) = 0.7 is the perturbation, since in the fixed point (73(0) = 0.0. As in the case g=0, we plotted a 
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set of phase diagrams and a graph of e versus SALI. In the SALI graphs, we break the vertical axis to show 
that the minimal SALI values are above the threshold of chaos, according to the expectations of the SALI 
method, that is, SALI < 1(T 5 . 




q 3 (0) = 0.7 and e varying as a)e = 0.10, b)e = 0.25, c) e = 0.50, d)e = 1.5, e) e = 2.0 and f ) e = 2.5. 

Figure 4 
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6 Concluding comments 

The comparison encompassing the integrable critical coupling "pure Chern-Simons" system, presented in 
the work of Rcf.[6], and our extended N=2-susy descent model, the latter deserving an integrability analysis 
procedure both for the general and critical coupling regimes, as reported here, suggests that the extra susy 
may be responsible [23] for the global non-intcgrability (in the strong Painlevc sense) situation found even 
in the C.S.-like regime. 

In studying chaos, two main possibilities have been checked. First, we have verified if non- vanishing values 
of g were able to drive the system (previously with initial conditions and parameters such that regularity 
was achieved for g = 0) into a chaotic regime. The second point we have tackled concerns the opposite effect, 
namely, if a chaotic configuration for g = may become regular whenever g becomes non-trivial. 

In the case of a non-critical coupling, as g increases from zero, a stable configuration for g = keeps its 
stability as g varies. In those cases, the difference is that, for larger values of g, the system becomes slightly 
more unstable, but its dynamics is still regular. Something similar happens for configurations that exhibit 
chaos for g = 0. In such cases, the system remains chaotic, but a little more unstable. For critical coupling, 
orbits that were chaotic for g = become now regular; this may indicate that the critical coupling plays the 
role of a stabiliser of our model. These results may be interpreted on the basis of two points: 

1. For the critical coupling regime, there occurs a partial decoupling between the variables <f> and M, and 
this reduces the non-linearity of the system. 

2. The quantity G(G = 1 — g 2 <t> 2 ) must be positive, with < G < 1. This must be so in order to ensure 
positivity of the energy, and the existence of a stable ground state. In Eqs. (16), (17) and (18), G 
accompanies all terms with time derivatives, and for large enough g or <j>, G becomes small, rendering 
the algebraic sectors of these equations dominant. This fact may have a stabilising consequence, 
implying that, in the case of a non-critical coupling, the dynamics for g ^ is not that different from 
the case with g = 0. In the critical coupling regime, the stabilising effect could be a combination of 
the two reasonings just presented. 

It is also noteworthy to point out that, for negative values of G, the phase space volume is no longer 
conserved, and, as a consequence, we do not have any longer a Hamiltonian system. For this reason, the 
results reported above make sense only if the initial conditions and the parameters ensure the positivity of 
G. Finally, as far as configurations with more regular behavior than the ones found in Ref.[6], as g increases, 
show up in this model, we wonder whether special physical conditions - extended supersymmctry in the 
original system - might be responsible for this stabilising effect. 
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